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Abstract 



A new formalism to calculate the transmission coefficient t of electrons from 
a material £ into the same material £ through a barrier B is presented. The 
barrier B is arbitrary and can be metallic, semiconducting or insulating. The 
important feature of this formalism is that it starts from the electronic struc- 
ture of a periodic C/B superlattice. The electronic stucture of such a superlat- 
tice is calculated selfconsistently within a Screened Korringa-Kohn-Rostoker 
method. The capability of the new method is demonstrated by means of a 
free electron model. First results for the transmission of Cu electrons through 
a Co layer are presented. 

I. INTRODUCTION 

If Bloch electrons encounter planar interfaces they may be reflected or transmitted. In 
this paper we develop a technique for the calculation of these reflection and transmission 
amplitudes within the Screened Korringa-Kohn-Rostoker (Screened KKR) approach, a tech- 
nique for ab initio electronic structure calculations. 

There are many uses for Bloch wave transmission and reflection amplitudes and proba- 
bilities. First of all they can be applied to describe transport jl],3 in films and multilayers 
including the electronic structure of the real material. Secondly, they can be used do develop 
a first principle based theory of electron tunneling. Finally, transmission and reflection co- 
efficients are the basic ingredients for the theoretical understanding of interlayer exchange 
coupling |3|]. 

Consider an arbitrary barrier B which is sandwiched by two leads of the material £ (see 
Fig. 1). For the propagating bulk-states in £ the transmission coefficient t through the 
barrier is defined as 

t = for B, = (1) 

where the Aq i and B\ are the complex amplitudes of the right- and left-moving Bloch waves. 
During the scattering process the energy E and the in-plane component kii of the wave vector 
{k, kii) are conserved. 

In the past different methods were established for the calculation of t, e.g. These 
methods usually model the structure of the two semi-infinite leads £ separated by a barrier 
B. Via the Green's function or the wave function of the system t can be computed. In 
this work we present a new method which approaches the problem differently: We use the 
band structure and wave function of a periodic superlattice B/£ to calculate t. This is 
strongly motivated by the development of highly efficient N-scaling superlattice KKR codes, 
e.g. M- 

A superlattice is a structure B/C which is periodically repeated in all three dimensions 
(see Fig. 2). Therefore the case of the single barrier B with leads £ is fully contained in this 
structure and the similarity with the Kronig-Penney model is obvious. Kronig and Penney 
showed in their simple 1-dimensional free electron model how the bandstructure of a periodic 
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system can be derived from the scattering properties of a single barrier. (For more details 
on the Kronig- Penney model see for instance 



The keynote of this paper is to derive the scattering properties of a single barrier from the 
electronic structure of a superlattice. We discuss the theory and numericalimplementation 
of the new method and present its first successful applications to free electrons and a Cu/Co 
(001) system. 



II. THEORY 

First we state our basic assumptions. Thereafter we discuss the transmission coefficient 
t of the single barrier and finally its connection to the superlattice wave function. 



A. Basic assumptions 

In Fig. 1 and Fig. 2 schematic pictures of a single barrier and a superlattice are presented. 
Note that x is always perpendicular to the barrier £>, the other two components are denoted 
with rii. Along rn the leads £ and the barrier B show the same 2-dimensional periodicity. 
The superlattice, however, is also periodic in x-direction. 

The single barrier case (see Fig. 1) consists of the barrier and two semi-infinite leads. In 
difference to the common understanding, we will include the perturbed interface regions of 
the leads into the barrier B. The leads C, however, consist of bulk-like potentials only. A 
superlattice (see Fig. 2) consists of an infinite number of barriers separated by interlayers. 
The barriers B contain again the perturbed interface regions. The interlayers consist of 
bulk-like potentials, the same as in the leads L of the single barrier case. 

Furthermore, we assume energy- and ^-conservation due to the elastic scattering and the 
in-plane symmetry of the system. Since kn is conserved the problem can be treated as quasi 
1-dimensional, but the reader should be cautious while keeping track of the components of 
the other two dimensions. 

We make two additional simplifications: First, the barrier is symmetric and, secondly, we 
restrict our considerations to one propagating bulk state $( fc ' k ll)(ic, rii) at a given energy and 
kn . That means, degeneration and band-crossing of bulk states in the leads are excluded. 



B. The single barrier 

The wavefunction in the leads £ is a superposition of the right- and left-moving bulk 
states, since the solution of the Schrodinger equation depends only on the local potential 
(see Fig. 1): 

tt(x,r,|) = ^e^e^ii-N^k^x^ii) + ^e-^e^H'n^k,, (rc, r,,) . (2) 

n = {0,1} denotes the two sides of the barrier. <3>( fc,k ii)(x, ry) = e lkx e lk ^ r ^Uk,k»(^, r \\) is the 
Bloch function of the bulk system. The decomposition of $^ fc,k n^(x, ry) into a phase factor 
and a periodic function is due to Bloch's theorem with u^,. (x, ry) = Uf.^,, (x — md, ry), where 
d is the bulk lattice parameter and m an integer. 
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The amplitudes A n and B n on both sides of the barrier are related to each other by a 
transfer matrix M which describes all the scattering processes inside the barrier: 

{$)-»{$)■ (3) 

The four complex variables of this M-matrix are not independent from each other. Due 
to time-reversal invariance, current conservation and the symmetry of the barrier, M can be 
parametrized with two independent real parameters a\ and (3\ in the following form [ |TT|j : 

M= f ai + ifk ifo \ (4) 

with (3l = a\ + Pl-\. 

Therefore the relation M _1 = M* is valid. 

Using the components of the transfer matrix the transmission coefficient t is by definition 
Eq. ©: 

* = — Kr ■ ( 5 ) 

Obviously, the problem is solved if the eigenstates ty(x,ru) of such a single barrier are 
known. 



C. The super lattice 

In contrast to the eigenfunctions of the single barrier (see Eq. @) the eigenfunctions of 
a superlattice ^(x, ry) are described by quantum numbers (k, ky). Here, k is the superlattice 
wave vector in x-direction reflecting the prolonged lattice structure with the lattice constant 
I, that is the length of the unit cell in x-direction. Furthermore, the in-plane component ky 
is the same as in the bulk system since in-plane symmetry is conserved. The translational 
symmetry of \I/( K ' k ii) in the superlattice direction is given by Bloch's theorem: 

^ k \\\x, ry) = (e~ iKl ) n ¥ K ' k \0(x + nl, ry) , (6) 

where the integer n ranges from minus to plus infinity. 

Let us now consider the superlattice as a periodic repetition of single barriers B with 
the interlayer £ in between. In the interlayer the potential is bulk-like and therefore we 
can expand the superlattice wave function (see Eq. (|])) as a superposition of right- and 
left-moving Bloch functions of the bulk system in the cell n = (a < x < a + s): 

A e jfc V k ii r ii Mfc , k|| (x,rn) + Boe- lfc VVii M _ fc ^(£,ry) (7) 
= [[e~ iKi y [A n e ik{x+nl) ^ r \\u KH (x + nl,r { \) + B n e" ik{x+nl) ^Wu.k^x + n/,ry)" . 

By using the relation Uk^Ax + nl, ry) = Uk^Ax, ry) we can rewrite this into 
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^tegikjr,, u r x rn)+B e ■■-< 



(e~ iKl ) n [A n e 



iknl ikx „ikiiri 



Uk,U K x , r ll 



-iknl —ikx „tkii r i 



M_ fc , k|l (x,r|i 



Since the right- and left-moving bulk states are orthogonal we find 



A „iknl 
J± n t 

td „—iknl 
J-> n t. 



(8) 
(9) 



(10) 



On the other hand we can connect the amplitudes in different cells with the transfer matrix 
M: 



A 
B 



A n 



(11) 



We can combine Eqs. ( |I0"D and (FTTT) to the eigenvalue equation 

= P 



inl 



A 
B 



A 
Bo 



(12) 



with 



P 



e ikl Q 

e~ ikl 



M' 1 



(cti - i(3i)e ikl -i(3 2 e ikl 
i(3 2 e- ikl {at + ip x )e~ ikl 



(13) 
(14) 



Since P is fully determined by the single eigenvalue e ird and the corresponding eigenvector 



Bo 



we obtain the following algebraic set for ot\ and fi\. 

(3 2 2 = a\ + ft - 1 

Ao = 

Bq ol\ sin kl — (3\ cos kl — sin kI 
cos kI = oti cos kl + /3i sin kl . 



(15) 



Thus, we are able to calculate t for a given energy and ky knowing k, k and ^ of the 
superlattice wave function \|>( K ' k ii) and corresponding bulk state $( fc ' k n). 



III. NUMERICAL ASPECTS 

In this section we explain how to extract k and ^ from \l>( K ' k ii^(x, ry) in the interlayer, 
which is very similar to a Fourier analysis of the ^^^{x, rii). For the electronic band 
structure calculations we used a Screened Korringa-Kohn-Rostoker (Screened KKR) method. 
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A. Brief introduction to the Screened KKR method 



The Screened KKR method (see J^-CT) was developed for the selfconsistent calculation 
of the electronic structure. It is based on Green's functions and the computation time 
scales linearly with the number of atoms in the unit cell. This allows the computation 
of systems with a relatively large number of atoms. For the atomic potentials the atomic 
sphere approximation (ASA) is used. The superlattice wave function is expanded in terms 
of spherical solutions: 

*«(R« + r M + r) e iKRn CZ( K )R?(r, E)Y L {r) , (16) 

L 

the wave function in the leads as well: 

$k (R n + ^ + r) = £ e *R*££(k)i¥(r, E )Y L {r) . (17) 

L 

R n is the lattice vector pointing to the n th unit cell, r M refers to the /j, th atom in the unit 
cell and r is the space coordinate inside the fj, atomic sphere. The Yl{t) are the spherical 
harmonics with the short hand notation L = (l,m), I is the angular momentum and m the 
magnetic quantum number. Rf (r, E) is the radial solution of the Schrodinger equation in 
the /j, th atomic sphere. The C^(n) and C^(k) are the expansion coefficients, k = (k, kn) and 
k =(k, k||) are the Bloch vectors in the superlattice and in the bulk system, respectively. 

Angular momenta up to / max = 3 were taken into account. 



B. Implementation 



As already explained, the superlattice wave function \I/( K ' k ii) in the interlayer C can be 
decomposed in right- and left-moving Bloch waves 



Yj/(«,k||) = y4 $( fc . k l|) + _g $(- fe ' k l|) . 



(18) 



Using the angular momentum expansions of Eqs. fll6|) and (|17|) we can expand Eq. (JlSf) in 
the th unit cell, which corresponds to R° = ( — Z, 0, 0): 



J2e- M C£(K,k {] )R?(r,E)Y L (r) 



E 



A e- m C^k,k ll )Rf(r,E)Y L (r) 
+B e ikl C£(-k, k|,)i?f (r, E)Y L (r) 



(19) 



Now we make use of the orthogonality relation of the Yi(r) and the following relations 
between the C£(±k, ky) due to the real potential of the leads and Bloch's theorem: 



C£(±Jfc,k||) = e i(±fc ' k ii )(r ^ ro) (7^(±A;,k||) 
^(-fc,k,|) = (-l)^(fc,k„) . 



(20) 
(21) 



So finally we conclude: 
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= e- w C^(«,k||) (22) 
-A e- ikl e i{ - k ^ ){r ^ ro) Cl(k^) 
-B e iH e i{ - k ' k ^- ro) C L {k, k||)(-l)' +m . 

We can exploit this equation at three different atomic positions \i in the interlayer and 



receive an algebraic set for k and the ratio It turns out that the C^(k, ky) cancel 



so 



no additional computation of the bulk eigenstates of the leads is necessary. All we have 
to compute are the eigenstates of the superlattice \&'" ,k ii'(:£, ry) for a given energy E and 
in-plane wave vector ky which gives us the C£(k, ky) and k. From here on we can calculate 



t via Eqs. ([151) and (p2]) . As a side effect we obtain the bandstructure k(E) of the lead 
material. This gives us the nice opportunity for an intermediate test of our method as we 
will see in the next section. 

Still there is another problem we are faced with: If k is complex the corresponding 
(x, ry) is evanescent which means it does not exist in the superlattice. In other words 
we are in a superlattice bandgap. But also in this energy range we want to calculate the 
corresponding k and t. One opportunity is the computation of the complex bandstructure 
and developing a similar formalism as the one described above. But usually the complex 
bandstructure is not available. Another nice idea is to vary the interlayer thickness s = 
(/ — 2a) (see Fig. 2): The superlattice bandstructure gets stretched or squeezed on the 
energy axis comparable to the eigenstates of a square well. In this way we can scan the 
whole energy range, albeit it makes more than one superlattice calculation necessary. We 
followed this second idea. 



IV. RESULTS 

The results we discuss in this section demonstrate that our new method is working 
successfully. We considered the scattering of free electrons through a rectangular barrier, 
which is analytically solvable. Furthermore we investigated the transmission of Cu electrons 
through 5 monolayers (ML) Co (001), which we compared with the results of Wildberger 
||12|| , who developed a different method for the computation of t. In addition to the \t\ (in- 
dependence we checked the dispersion relation k(E) of the lead material. 



A. Free electrons through a rectangular barrier 

Outside the barrier we modelled an exact free electron gas with constant potential V = 0. 
The barrier itself consists of five layers with muffin tin (MT) spheres of a constant potential 
inside the sphere. The potential height is V = 0.125Ryd and they are arranged in an fee 
lattice with a lattice constant of d = 6. 76a.it. . This is only an approximate rectangular 
barrier, but the best one possible with ASA potentials. 

In Fig. 3 the dependence of \t\ on (E — kjj) for ky = |0 is plotted. The analytical 
curve was fitted to the numerical data with a potential height of the rectangular barrier 
V = 0.76 • 0.125Ryd. This potential height is given by the dashed line. The factor 0.76 is 
very close to 0.74 which is the effective space filling of MT spheres on an fee lattice. The 
remaining small difference of 0.02 and also the subtle deviations in the \t\ ( independence can 
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be understood very well in terms of the different barrier shapes. The two numerical curves 
with the different interlayer thicknesses (circle and plus) coincide, which is a brilliant test for 
the validity of the method itself. The discontinuity of the curves is due to the superlattice 
bandgaps. Obviously these gaps can be filled with different interlayer thicknesses s. To fill 
all gaps more calculations with different interlayer thicknesses s would be necessary. This 
is not done here since we just want to demonstrate the capability of our new method. For 
a short interpretation we notice that the electrons have a low transmission coefficient below 
the barrier since they have to tunnel through. Above the barrier height they approach t — 1 
very rapidly, the oscillations are due to interference effects on the barrier interfaces. 

In Fig. 4 we show the correspondence between the calculated dependence of k on (E — k|) 
for different ky (circle, cross and triangle) and the expected free electron dispersion (solid 
line), and obtain a very good agreement. 

B. Cu electrons through 5 ML of Co (001) 

We arranged the Co and Cu on a fee lattice with a lattice constant of 6. 76a. it.. The value 
is larger than that of magnetic fee Co and smaller than that of fee Cu. The Co layer consists 
of 5 ML. The potentials of the Co layer and three Cu layers at both sides of the barrier 
are taken from a self-consistent C05CU7 supercell calculation and form the barrier B. The 
potential in the middle of the Cu layer is a bulk-like Cu ASA potential. This approximation 
has no influence on the results and has the advantage that we can construct easily any Cu 
interlayer thickness s without doing an extra selfconsistent calculation. 

Figure 5b and 5d show the \t\ (E')-dependence for ki = in the majority and minority 
channels. Since Co is a ferromagnetic material the superlattice bandstructure of Cu/Co has 
no more spin-degeneracy and we obtain spin-splited majority and minority bands. As for 
the free electron results the curves for different interlayer thicknesses s match, which is not 
shown explicitly here. Large parts of the \t\ (-E)-curve can be justified qualitatively by the 
band structure (BS) of Co and Cu, which are plotted in Fig. 5a (Co majority BS), Fig. 
5c (Cu BS) and Fig. 5e (Co minority BS). First we discuss the \t\ ( Independence in the 
majority channel (see Fig. 5b): 

Since the Co and Cu states have the same symmetry character in the energy region I 
the Cu electrons can propagate through the barrier very well, below the bottom of the Co 
band they have to tunnel. Therefore the behavior in section I is very similar to the free 
electron result as discussed in the previous subsection. In region II there are no propagating 
Cu bulk states and consequently there is no t. In region III we see that t is almost free 
electron like again, the rapid breakdown at about 0.45 Ryd corresponds the change of the 
symmetry character in the Co bands. Region IV has to be excluded from our discussion 
since the Cu band structure is degenerated there and hence our method is not applicable. 
In region V, including the Fermi level at 0.68Ryd, we obtain a £ of almost 1 which is due to 
the corresponding band symmetries in Co and Cu. 

In the minority channel (see Fig. 5d) we see two major differences to the majority channel 
in region III and V: In region III t is almost zero since the Co minority band has a gap 
there and so no band matching is possible between the Cu and Co states. In region V we 
see a low t below 0.64 Ryd. This is due to the upward shifted Co bands in comparison to the 
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majority channel and therefore no band matching between states with the same symmetry 
character occurs. 

Furthermore, we compare our results in the energy region V (see Fig. 6a) with the results 
obtained with a different method by Wildberger |12| (see Fig. 6b), who used a semi-infinite 
structure with two Co barriers at a certain distance (which is comparable to our interlayer 
thickness s) for the calculation of the reflection coefficient |r|. The reflection coefficient \r\ is 

related to \t\ by \r\ = v/l — \t\ 2 . In Fig. 6b we see their \r\ ( independence for the majority- 
and minority-channel (dashed and solid line) for kii = around the Fermi energy. The 
agreement with our own results shown in Fig. 6a is very good. Slight differences are related 
to the different lattice parameters which were used. In addition, Wildberger evaluates the 
Green's function of the system for energies with a small imaginary part, which causes a 
broadening corresponding to a finite temperature. The agreement between our calculation 
and the calculation of Wildberger is a confirmation of the capability and the validity of our 
new method. 

In Fig. 7 we see the band structure of fee Cu along the line T — X. The solid line is 
taken from a pure Cu band structure calculation, the crosses denote the band structure k(E) 
extracted from the superlattice calculation. The region between the dashed lines corresponds 
to region IV in Fig. 5 and has to be excluded for the reasons mentioned above. Again, the 
correspondence between both curves is perfect. 

V. SUMMARY 

We presented a new method to calculate the transmission coefficient t of bulk states 
through a single barrier from a superlattice band structure. Its validity is fully confirmed 
by the successful test results. At this state the method is restricted to symmetric barriers 
and to only one propagating bulk state at a given energy and in plane wave vector. The 
generalization to the case of mor than one propagating bulk state and unsymmetric barriers 
is under development. 
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FIG. 1. Schematic drawing of the single barrier case with a barrier B (including the perturbed 
interface regions) and leads C. Furthermore, right- and left-moving eigenstates <J> are shown. 
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FIG. 2. Schematic drawing of the superlattice with barriers B (including the perturbed inter- 
face regions) and interlayer C. The superlattice wave function in the interlayers is a superposition 
of right- and left-moving bulk eigenstates $>. 
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FIG. 3. The dependence of \t\ on (E - kg) for ky = of free electrons through a 

rectangular barrier with a thickness of 5 ML and an averaged height of 0.092 Ryd denoted by the 
dashed line: Analytical results (solid line) in comparison with the results obtained with our new 
method for two different interlayer thicknesses (cross: 11 ML, circle: 17 ML). 
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FIG. 4. The dispersion relation of free electrons for different kii (circle: k x = k y = 0, cross: 
k x = 0.5, k y = 0; triangle: k x = k y = 0.5, all in relative units of the Brillouin zone) in comparison 
with the analytical result (solid line). 
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FIG. 5. Band structure of the majority (a) and minority (e) electrons of Co, band structure of 
Cu (c), furthermore the \t\ (-E')-dependence for Cu electrons through the majority (b) and minority 
(d) potential barrier of 5 ML Co; all along the line T — X, that is ky = 0. 
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FIG. 6. The dependence of |r| on (E — Ep) with our new method (a) for the majority (dashed 
line) and minority (solid line) Cu electrons through 5 ML of Co(100) at kii = in comparisson 
with Wildberger (b) (see Fig. 8.5, p. 152 in G2|); the axes are the same for both figures. 
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FIG. 7. Band structure of fee Cu along FX from a Cu bulk calculation (solid line) and 
indirectly from the superlattice Cu/Co (001) calculation (cross). 
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